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Abstract-- A possible solution of the well known paradox of chiral molecules is based on the 
idea of spontaneous symmetry breaking. At low pressure the molecules are delocalized between 
the two minima of a given molecular potential while at higher pressure they become localized 
in one minimum due to the intermolecular dipole-dipole interactions. Evidence for such a phase 
transition is provided by measurements of the inversion spectrum of ammonia and deuterated 
ammonia at different pressures. In particular, at pressure greater than a critical value no inversion 
line is observed. These data are well accounted for by a model previously developed and recently 
extended to mixtures. In the present paper, we discuss the variation of the critical pressure in 
binary mixtures as a function of the fractions of the constituents. 


1. INTRODUCTION 

According to quantum mechanics chiral molecules, that is, molecules which are not superimposable 
on their mirror image, should not exist as stable stationary states. Consider ammonia NH 3 . The 
two possible positions of the N atom with respect to the plane of the H atoms are separated by a 
potential barrier and can be connected via tunneling. This gives rise to stationary wave functions 
delocalized over the two minima of the potential and of definite parity. In particular, the ground 
state is expected to be even under parity. Tunneling induces a doublet structure of the energy 
levels. 

On the other hand, the existence of chiral molecules can be interpreted as a phase transition. 
In fact, isolated molecules do not exist in nature and the effect of the environing molecules must 
be taken into account to explain phenomena characterized by instabilities. This interpretation 
underlies a simple mean-field model developed in [1] to describe the transition of a gas of NH 3 
molecules from a nonpolar phase to a polar one through a localization phenomenon which gives 
rise to the appearance of an electric dipole moment. Even if ammonia molecules are only pre-chiral 
[ 2 ], the mechanism, as emphasized in [T], provides the key to understand the origin of chirality. 

A quantitative discussion of the collective effects induced by coupling a molecule to the envi¬ 
ronment constituted by the other molecules of the gas was made in [3]. In this work it was shown 
that, due to the instability of tunneling under weak perturbations, the order of magnitude of the 
molecular dipole-dipole interaction may account for localized ground states. This suggested that a 
transition to localized states should happen when the interaction among the molecules is increased. 

Evidence for such a transition was provided by measurements of the dependence of the doublet 
frequency under increasing pressure: the frequency vanishes for a critical pressure P CI different for 
NH 3 and ND 3 . The measurements were taken at the end of the 1940s and beginning of the 1950s 
El El E] but, as far as we know, no quantitative universally accepted theoretical explanation exists 
in spite of many attempts. The model jlj gives a satisfactory account of the empirical results. A 
remarkable feature of the model is that there are no free parameters. In particular, it describes 
quantitatively the shift to zero-frequency of the inversion line of NH 3 and ND 3 on increasing the 
pressure. 

Recently, we extended the model [lj to gas mixtures [7|. This case may be of interest, among 
other things, for the interpretation of the astronomical data such as those from Galileo spacecraft 
| 8 j which measured the absorption spectrum of NH 3 in the Jovian atmosphere. Formulas for the 
critical pressure of a general mixture have been provided. 

In the present paper, we investigate the behavior of the critical pressure in binary mixtures. We 
show that the critical pressure of a chiral species can be increased or decreased by several orders 
of magnitude by mixing it with a proper fraction of a proper species, chiral or non polar. 
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2. CHIRAL GAS 

We model a gas of all equal chiral molecules as a set of two-level quantum systems, that mimic 
the inversion degree of freedom of an isolated molecule, mutually interacting via the dipole-dipole 
electric force. At moderate density, we approximate the behavior of the N 3> 1 molecules of the 
gas with the mean-field Hamiltonian 

AE -x 

~T a - ' 

where a x and <r z are the Pauli matrices and ip is the Pauli spinor representing the mean-field 
molecular state with normalization (ip, ip) = N. The scalar product between two Pauli spinors is 
defined in terms of their two components in the standard way. The parameter A E is the inversion 
energy-splitting measured by spectoscopic methods in the rarefied gas. The parameter G accounts 
for the effective dipole interaction energy of a single molecule with the rest of the gas. It can be 
estimated in two different but equivalent ways [7j. 

The first way to estimate G is based on the so called Keesom energy, namely, G is identified with 
the effective dipole-dipole interaction obtained after averaging over all possible molecular distances 
and all possible dipole orientations. These averages are calculated assuming that, concerning the 
translational, vibrational, and rotational degrees of freedom, the N molecules behave as an ideal gas 
at thermal equilibrium at temperature T. This assumption relies on a sharp separation (decoupling) 
between these degrees of freedom and the inversion motion. The result is 
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where P is the pressure of the gas, p the electric-dipole moment of the molecules and d the mini¬ 
mal distance between two molecules, namely, the so called molecular collision diameter. At fixed 
temperature, the effective interaction constant G increases linearly with the gas pressure P. 

The second way to estimate G is based on the reaction field mechanism [9]. Let us consider a 
spherical cavity of radius a in a homogeneous dielectric medium characterized by a relative dielectric 
constant s r . An electric dipole p placed at the center of the cavity polarizes the dielectric medium 
inducing inside the sphere a reaction held R proportional to p. As a result, the dipole acquires an 
energy 
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Since e r ~ 1, we can approximate the first fraction in Eq. <§ by the Clausius-Mossotti relation 
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where a is the molecular polarizability and oP = p 2 /(?>£$kgT) is the Debye (orientation) polariz¬ 
ability. Observing that for a chiral gas a D 3> a (for instance, in the case of NH 3 we have a ~ 2 A 
whereas a D ~ 217 A 3 at T = 300 K), we get 
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Microscopic arguments mm show that the radius of the spherical cavity a is not arbitrary but 
must be identified with the the minimum distance between two interacting molecules, namely, the 
molecular collision diameter d introduced in Eq. S. We thus have $ = —G. 

The state ip collectively describing the inversion degree of freedom of the gas of N chiral molecules 
is determined as the minimal-energy stationary state of the Hamiltionan of Eq. (|T|. This corre¬ 
sponds to find the lowest-energy eigenstate of the nonlinear eigenvalue problem h(yj)ip = A ip, with 
the constraint (ip, ip) = N. We refer the reader to [l| Uj for the mathematical details, here we just 
state the main results. 

There exists a critical value of the interaction strength, G cr = AE/2, such that for G < G cr 
the mean-held eigenstate with minimal energy is ip = \/N<p+, where </?+ is the eigenstate of a x 
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with eigenvalue +1, i.e., the molecules are delocalized. For G > G cr there are two degenerate 
solutions of minimal energy which can be termed chiral states, in the sense they are transformed 
into each other by the parity operator a x . For G G cr , these solutions become the localized states 
y/Ni^L, VN<Pr, where <pl, Pr are the eigenstates of a z . The energy associated with the state i/j is 
a continuous function of G with a discontinuous derivative at G = G a -. We thus have a quantum 
phase transition between a delocalized (or achiral, or nonpolar) phase and a localized (or chiral, or 
polar) phase. In view of the dependence of G on P, we can define a critical pressure at which the 
phase transition takes place 
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Note that the value of P CT is completely determined in terms of the microscopic parameters A E, p 
and d and the temperature T. 

In [1] we have also shown that in the delocalized phase the inversion angular frequency of 
the interacting molecules depends on the pressure as hu(P) = AEy/l — P/P CT . This formula is 
interesting as it expresses the ratio of two microscopic quantities, hu and A E, as a universal function 
of the ratio of the macroscopic variables P and P cr . Furthermore, it is in very good agreement with 
some spectroscopic data showing the shift to zero frequency of the inversion line of NH 3 or ND 3 at 
increasing pressures 01 El El. 


3. BINARY MIXTURES 

Consider a gas mixture of two species labeled 1 and 2. In this case, the Clausius-Mossotti relation 
reads as 

£ _|_ 2 = 3 (d 1 ( ai + “?) + P' 2 (“2 + a 2 ) )) 1 (7) 

where the Debye polarization a? = p 2 / (3£ok B T) is given in terms of the molecular electric-dip ole 
moment m of the species i = 1,2. According to the reaction field arguments discussed above, a 
chiral molecule, let us say of species i, having dipole moment p n . acquires, due to the interaction 
with all the other molecules of the mixture, the energy 

1 u 2 

S’i = —- (pi (aq + a?) + p 2 (a2 + a?)) ^7-^3 > (8) 

where d% is the molecular collision diameter of the i-th species. We now specialize the discussion in 
the following two cases. 

3.1. Mixtures of a chiral gas with a non polar gas 

For a mixture of chiral and non polar molecules, the mean-field molecular state of the chiral species, 
assumed as species 1, is determined similarly to the case of a single chiral gas. The only degree of 
freedom of the non polar molecules is the deformation which, in turn, is proportional to the electric 
dipole moment of the chiral molecules. We may thus describe the mixture by a single mean-field 
molecular state, ipi, normalized to the number of molecules of the species 1, (t/x= N±. As 
before, we assume that this state is determined as the lowest-energy eigenstate of the eigenvalue 
problem associated with the mean-field Hamiltonian 
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In this Hamiltonian —G\ represents the effective dipole interaction energy of a single chiral molecule 
with all the other molecules, chiral and non polar, of the mixture. Thus we can identify —G 1 = S’i, 
where S\ is given by Eq. (J 8 | with i = 1, ot\ <C af, and a J? = 0, namely, 

G\ = (711 .Pi + 712 P 2 ), (10) 
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Figure 1: Critical pressure of the localization phase transition in a binary mixture of NH 3 as a function of 
the fraction of the second constituent chosen as ND3 (bottom-left axes) or He (top-right axes). 


As usual we used the ideal gas relations pi = Pi/ksT and p 2 = Th/fesT, where P\ and P 2 are the 
partial pressures of the two species. 

The analysis of the nonlinear eigenvalue problem h{^i)ipi = Ai?/q, with the constraint {i/j 1 , Vh) = 
Ni, is identical to the case of a single chiral gas. We have a localization phase transition when 
G\ = AEi/2. The transition can be considered as a function of the total pressure P = P\ + P 2 of 
the mixture and of the fractions of the two species x\ = P\/P and X 2 = P-ijP ■ In this case, instead 
of a unique critical pressure, we have a critical line parametrized by x\ or X 2 = 1 — xi, e.g., 
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In Fig. [lj we show the variation of the critical pressure in a NH^-He mixture as a function of the He 
fraction. The value of P cr increases from the critical pressure of pure NH 3 to a maximum reached 
for a vanishing NH 3 fraction. The value of this maximum depends on the nature of the non polar 
species, it is greater the smaller is the molecular polarizability 07 - 

3.2. Mixtures of two chiral gases 

For a mixture of two chiral gases, we describe the inversion degrees of freedom of the two species 
by mean-field molecular states V , i,V ; 2 normalized to the number of molecules of the corresponding 
species. These states are obtained as the lowest-energy eigenstates of the eigenvalue problem 
associated with the coupled mean-field Hamiltonians 
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Pauli operators now have a label i = 1,2 relative to the species they refer to. Each term —Gij 
represents the effective dipole interaction energy of a single molecule of species i with all the other 
molecules of species j. By matching —Gn — Ga = S’i, where Si is given by Eq. ([ 8 | with ai <S a ^ 
and 0:2 <C a®, we get 
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where P is the total pressure of the mixture and x\. X 2 the fractions of the components. 

The solution of the coupled nonlinear eigenvalue problem hi(ipitp 2 )ipi = Ai^i and ^(V’l^ 2)^2 = 
X 21 P 2 , with the constraints (ip\,ip\) = N\ and (ip 2 , 1 ( 2 ) = -/V 2 , is discussed in [7]. As in the case 
of a single chiral gas, the mixture undergoes a localization phase transition at a critical pressure 
P cr . For 0 < P < -P cr , the lowest-energy molecular state of the mixture corresponds to molecules of 
both species in a delocalized symmetric configuration. For P > P cr , new minimal-energy molecular 
states appear with twofold degeneracy. These states correspond to molecules of both species in a 
chiral configuration of type L or R. We have a particularly simple formula for P cr , 
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the inverse critical pressure of the mixture is the fraction-weighted average of the inverse critical 
pressures of its components. In Fig. [TJ we show the variation of the critical pressure in a NH 3 -ND 3 
mixture as a function of the ND 3 fraction. The value of P cr ranges from the critical pressure of pure 
NH 3 to that of pure ND 3 , namely, from 1.69 atm to 0.11 atm. By changing ND 3 with, for instance, 
D 2 S 2 , the minimal value of P cr can be extended down to 4.3 x 10 ~ 9 atm, the critical pressure of 
deuterated disulfane. 


4. CONCLUSION 

Our approach to the existence of chiral molecules is based on ideas of equilibrium statistical mechan¬ 
ics. One may be surprised by the presence of a quantum phase transition at room temperatures. 
We emphasize that the transition takes place only in the inversion degrees of freedom. The dynam¬ 
ics of these degrees of freedom is affected by temperature only through the values of the coupling 
constants. 

We have shown that with the addition of a proper fraction of a second species, non polar 
or chiral, the critical pressure of an ammonia mixture can vary in a range of several orders of 
magnitude. As a consequence, the inversion line of ammonia, as well as, possibly, that of the 
second chiral constituent, should undergo a frequency shift rather different from that measured for 
pure gases, see [7]. An experimental verification of these predictions is well within the reach of 
present technology and would represent a critical test of our theory. 
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